On nonclassical symmetries of generalized Huxley 
equations 

N.M. Ivanova^^ and C. Sophocleous^ 

^ Institute of Mathematics of NAS of Ukraine, 3 Tereshchenkivska Str., 01601 Kyiv, Ukraine 
e-mail: ivanova@imath.kiev.ua 

^ Department of Mathematics and Statistics, University of Cyprus, CY 1678 Nicosia, Cyprus 
e-mail: christod@ucy.ac.cy 

Nonclassical symmetries of a class of generalized Huxley equations of form ut — u^x + 
k{x)u'^{\ — u) are found. More precisely, for the class under consideration we completely 
classify reduction operators with r = 1 and give a wide number of examples of equations 
admitting reduction operators with r = 0. 

1 Introduction 

Consider reaction-diffusion equation of form 

ut = Uxx + k{x)u^{l-u), (1) 

where k{x) ^ 0. This equation models many phenomena that occur in different areas of math- 
ematical physics and biology. In particular, it can be used to describe the spread of a recessive 
advantageous allele through a population in which there are only two possible alleles at the 
locus in question. Equation ([T|) is interesting also in the area of nerve axon potentials [18j. Case 
k = const is the famous Huxley equation. For more details about application see l3^L4j and 
references therein. 

Bluman and Cole [2] introduced a new method for finding group-invariant (called also simi- 
larity) solutions of partial differential equations. The method was called by the authors "non- 
classical" to emphasize the difference between it and the "classical" Lie reduction method de- 
scribed, e.g., in jl4pi5). A precise and rigorous definition of nonclassical invariance was firstly 
formulated in |10j as "a generalization of the Lie definition of invariance" (see also |21j). Later 
operators satisfying the nonclassical invariance criterion were also called, by different authors, 
nonclassical symmetries, conditional symmetries, Q-conditional symmetries and reduction oper- 
ators [8l[9l[13]. The necessary definitions, including ones of equivalence of reduction operators, 
and relevant statements on this subject are collected in [TTlfTU]. 

Bradshaw-Hajek at al |3j;4J started studying class ([T|) from the symmetry point of view. More 
precisely, they found some cases of equations ([T]) admitting Lie and/or nonclassical symmetries. 
Complete classification of Lie symmetries of class ([TJ is performed in |12] . Conditional symme- 
tries of Huxley and Burgers-Haxley equations having nontrivial intersection with class ([T]) are 
investigated in [HI5H71ITT] . The present paper is a step towards to the complete classification of 
nonclassical symmetries of class ([T]). 

Reduction operators of equations ([I]) have the general form Q = rdt + i,dx + r]du, where r, ^ 
and rj are functions of t, x and u, and (r, ^) / (0,0). We consider two cases: 

1. r ^ 0. Without loss of generality r = 1. 

2. T = 0, ^ / 0. Without loss of generality ^ = 1. 

We present a complete classification for the case 1, while for the case 2 we found several examples. 



2 Equivalence transformations and Lie symmetries 



Since classification of nonclassical symmetries is impossible without detailed knowledge of Lie 
invariance properties, we review [12] the equivalence group and results of the group classification 
of class ([I]) . The complete equivalence group G~ of class ([T]) contains only scaling and translation 
transformations of independent variables t and x. More precisely it consists of transformations 

t = e\t + 62-, x = eix + e3, u = u, k = e^'^k, 

where Si, i = 1, 2, 3 are arbitrary constants, ei ^ 0. 

Theorem 1. There exists three -inequivalent cases of equations from class ([1]) admitting 
nontrivial Lie invariance algebras (the values of k are given together with the corresponding 
maximal Lie invariance algebras, c = constj (see 112]): 

1 : Vfc, {dt); 

2: k = c, {dt, d^); 

3: A; = CT~^ {dt, 2tdt + xd^) . 

In the following section we search for nonclassical symmetries which are not equivalent to the 
above Lie symmetries. 

3 Nonclassical symmetries 

Now let us recall the definition of nonclassical symmetry (or conditional symmetry, or reduction 
operator). 

Reduction operators (nonclassical symmetries, Q-conditional symmetries) of a differential 
equation C of form L{t,x,U(^r)) = have the general form 

Q = Tdt + + rjdu, 

where r, ^ and r] are functions of t, x and n, and (r, ^) ^ (0,0). Here denotes the set of all 
the derivatives of the function u with respect to t and x of order not greater than r, including 
u as the derivative of order zero. 

The first-order differential function Q[u] := rj{t,x,u) — T{t,x,u)ut — ^{t,x,u)ux is called the 
characteristic of the operator Q. The characteristic PDE Q[u] = is called also the invariant 
surface condition. Denote the manifold defined by the set of all the differential consequences of 
the characteristic equation Q[u] =0 in the jet space 

Definition 1. The differential equation C of form L{t, x,U(^j.)) = is called conditionally (non- 
classicaly) invariant with respect to the operator Q if the relation Q(r)L(i, x, ti(r))|^p|g(r)= 
holds, which is called the conditional invariance criterion. Then Q is called an operator of 
conditional symmetry (or Q-conditional symmetry, nonclassical symmetry, reduction operator 
etc) of the equation £. 

In Definition [1] the symbol Q(^j.) stands for the standard r-th prolongation of the operator Q 

[ails]. 

The classical (Lie) symmetries are, in fact, partial cases of nonclassical symmetries. Therefore, 
below we solve the problem on finding only pure nonclassical symmetries which are not equivalent 
to classical ones. Moreover, our approach is based on application of the notion of equivalence of 
nonclassical symmetries with respect to a transformation group (see, e.g., [IT])- more details, 
necessary definitions and properties of nonclassical symmetries we refer the reader to [16 t [T7 1 [T9 l 
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Since ([T|) is an evolution equation, there exist two principally different cases of finding Q: 
1. T / and 2. r = 0. 

First we consider the case with r 7^ 0. Here without loss of generality we can assume that 
r = 1. The results are summarized in the following theorem. 

Theorem 2. All possible cases of equations ([1]) admitting nonclassical symmetries with r = 1 
are exhausted by the following ones: 

1. k = c tan^ x: Q = dt — cot xdx , 

2. k = ctanh^ x: Q = dt — coth xdx, 

3. k = ccoth^ x: Q = dt — tanhxdx, 

4. k = cx^: Q = dt - ^dx, 

5. A; = 4 (c> 0); Q = dt ± f (3m - l)dx - ^u^{u - l)du 

6. k = 2x-2; Q = dt + l{u- l)dx - ^u{u - Ifdu, 
where c is an arbitrary constant. 

Proof. We search for reduction operator (operator of nonclassical (Q-conditional) symmetry) 
in form Q = dt + ^{t, x,u)dx + r](t, x,u)du. Then the system of determining equations for the 
coefficients of operator Q has the form 

Cuu — 0, 

2^^^ - 2ri^u - ^KuU^ + ^KuU^ + 2r]xu - ixx + = 0, 

-kr]uU^{l -u) + 2/c^^M^(l -u)+ r]xx - 2^xV 

= Vt — kxiu^{l — u) — 2kriu + 3kr]u^ . 
From the first equation we obtain immediately that 

^ = (f){t,x)u + ip{t,x). 
Substituting it to the second equation we derive 

r] = -^(t>^u^ - (jiilJU^ + (t)xU^ + A{t, x)u + B{t, x). 

Then, splitting the rest of determining equations with respect to different powers of u implies 
the following system of equations for coefficients (/>, ij:, A and B. 
2 

-cj)^ - 3k(f) = 0, -4:(j)4>x + + 3k(f) = 0, 
3 

-2<Px^P + <Pt- 2<pipx - 2(pA + Scpxx = 0, 

2V'Vx - 20B + 2Ax - ^pxx + V't = 0, 

1 2 1 2 
--^4> 4>x + -k(p + k(f)tl! = 3k(j)x + kx(t), 

2 2 8 

-(/.Vx - - -^(kl - 2A:V'x - ^kA + 2#^.V' 

2 

= --^(p<Pt - '2.k<l)x - kx(f) + kx-ip, 

-2<l)xA + (pxxx + 2(/'V'V'x + kA- (pxxip - Hx^Jx - (pi'xx + "^kipx 
= 4>tx - - (P'lpt - kxip + 3kB, 
Axx - 2iPxA = At + 2(t)xB - 2kB, 

-2ijxB + Bxx = Bt. (2) 
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Now from the first equation of ([2]) it is obvious that either (i) = or (ii) (pt = 0, k = |(/)^. 
Consider separately these two possibihties. 



Case (i). (p{t,x) = 0. System ([2]) is read now hke 

2A^ - + 2'0Vx + V't = 0, 
-2/cV'x - '2kA - ipk^ = 0, 
2/cVx + kA + k^Tp - 3kB = 0, 
A^^ - 2AiP^ + 2kB -At = Q, 
B^^ - 2tP^B -Bt = 0. 

From the second and third equations we deduce that A = —3B. After substituting this to the 
previous system we obtain 

2ipip^ - 6Bx - ipxx + V"* = 0, 
kxi^ + 2k'4)x - 6kB = 0, 
GipxB - 3Bxx + 3Bt + 2k B = 0, 
Bxx - 2SVx -Bt = 0. 

It follows from the last two equations that B = 0. Then the rest of the determining equations 
is read like 

27pipx - tpxx + = 0, kxtp + 2kil)x = 0. 
General solution of this system is 

k = c, = const, 

c , ax + b 

^ = 7 TTTI ' ^ = ^+ — [ ' 

[ax + oj^ 2ta + m 

k = ip' = tp'^ + a. 

tp"^ 

Nonclassical symmetry operator obtained from the first two branches of the solution of the 
above system are equivalent to the usual Lie symmetry. The third branch (up to equivalence 
transformations of scaling and translations of x) gives cases 1-4 of the theorem. 

Case (ii). pt = 0, k = |(/)^. Substituting this to system ([2]) we obtain easily that ipt = At = 
Bf = 0. Then, the rest of the system ([2]) has the form 

-4Px + 2PiP + ^p^ = 0, 

-2Pxi^ - 2'pxP - 2pA + 3pxx = 0, 

2VVx - 2PB + 2Ax - tpxx = 0, 

2 2 8 4 14 8 

-P'^lpx - -4>4>xx - -Pl - -P'^A + y(/>V"Ax = --(p'^Px, 

4 2 2 

Pxxx + -P'^lpx - -4>'^B + -P>'^A - (plpxx - <Pxx1p - Hxlpx - 2p)xA + 2p'lp-px 

y o y 

4 

= --^PPxtp - P>'Pt, 

Axx - 2^PxA = -^(P^B + 2PxB + At, 
Bxx = 2ijxB + Bt. 
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It follows then that 

2(202 _ 9^^) ' 3</. ' 6(A ■ 

Substituting this values to the above system we obtain a system of 4 differential equations for 
one function <p only that has first order differential consequence of form 3(j)'^ + (^'^(px = 0. It is 
not difficult to show that general solution of this constraint (p = , (f> = c satisfies the whole 
system for (p. These two values of (p (taken up to equivalence transformations) give respectively 
cases 6 and 5 of the theorem. □ 

Note 1. Cases 1, 2 and 4 with c > were known in [Sjll], constant coefficient case 5 with c = 2 
can be found in, e.g., [llj . while 1,2 and 4 with c < 0, 3 and 6 are new. 

Now we turn into the case 2. That is, we consider nonclassical symmetry operator with r = 0. 
Without loss of generality we assume it to be of the form 

Q = dx + r]{t,x,u)du- 

Any operator nonclassical symmetry of the above form satisfies the following equation 
-r]xx - '^VVxu - ffriuu + kriuU^ - kijuU^ + ??t + kxU^ - kxV? - 2kr]u + 3ki]u^ = 0. (3) 

The corresponding invariant surface condition is Ux = i]. Eliminating Ux and Uxx^ equation ([T]) 
reads 

ut = ri'qu + 'nx + k{x)v?{l-u). (4) 

Using a solution of (l3|), we can find u by integrating invariance surface condition and then 
substituting in (jl]) to derive a solution of 

Now as it was shown in [16p20] for more general case of (l+n)-dimensional evolution equations, 
integration of equation ([3]) is, in some sense, equivalent to integration of the initial equation ([T]). 
However, since it contains bigger number of unknown variables it is possible to construct certain 
partial solutions. Thus, for example, we have succeeded to find all (C-inequivalent) partial 
solutions of equation ([3]) of the form 

n 

r]{x,t,u) = ^ (pp{x,t)vP, 

p=—m 

where m and n are positive integers and (pp{x,t) unknown functions. We find the following 
results: 

1. k = 2B^, T] = B{x)v?' — tanxn, where 

-4BB' + 4B' tan x-B" + 2B + 2B'^ tan x = 0. 
For example, this equation has a solution of form B = tanx. 

2. k = 2B^, 7] = B{x)v?' +tanhxu, where 

ABB' + AB' tanh x + B" + 2B + 2B'^ tanh x = 0. 
In particular, this equation has a solution of form B = — tanhx. 
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3. k = 2B^, T] = B{x)u'^ + cothxti, where 

ABB' + 4B' coth x + B" + 2B + 2B^ coth x = 0. 
In particular, this equation has a solution of form B = — cothx. 
A. k = 2^2, r] = B{x)u^ + |, where 

AxBB' + 45' + xB" + 2^2 = 0. 

A solution of this equation is B = — ^. 

5. A: = ^,r/ = i(n2-l). 

6. k = ^,r]=^u'^. 

7. k = 2 tan^ 2x, r/ = — n-^ tan 2x. 

8. A; = 2 tanh^ 2x, r] = tanh 2x. 

More detailed investigation of conditional symmetries and construction of associate similarity 
solutions of equations from class ([T]) will be the subject of a forthcoming paper. 
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